Hubbard exciton revealed by time-domain optical spectroscopy 
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We use broadband ultra-fast pump-probe spectroscopy in the visible range to study the lowest 
excitations across the Mott-Hubbard gap in the orbitally ordered insulator YVO3. Separating 
thermal and non-thermal contributions to the optical transients, we show that the total spectral 
weight of the two lowest peaks is conserved, demonstrating that both excitations correspond to the 
same multiplet. The pump-induced transfer of spectral weight between the two peaks reveals that 
the low-energy one is a Hubbard exciton, i.e. a resonance or bound state between a doublon and a 
holon. Finally, we speculate that the pump-driven spin-disorder can be used to quantify the kinetic 
energy gain of the excitons in the ferromagnetic phase. 



The physical properties of materials characterized by 
strong electron-electron interactions are determined by 
the competitive minimization of the potential energy and 
kinetic energy of the electrons. While the potential en- 
ergy is lowest for localized electrons, the tendency to a 
metallic state increases it at the expense of kinetic en- 
ergy. The fine tuning of the system's parameter results 
in the exotic ordering phenomena characterizing transi- 
tion metal oxides' (TMOs) phase diagrams[9]. In a Mott- 
Hubbard insulator the lowest electronic excitation across 
the gap creates, in the most simple case, an empty site 
(holon in the lower Hubbard band) and a doubly occu- 
pied site (doublon in the upper Hubbard band) [TO] . In 
the single-band Hubbard model, the energy of this tran- 
sition is solely determined by the on-site Coulomb re- 
pulsion between electrons leading to an effective energy 
cost U. Typically, holon and doublon are not bound 
to each other, but more composite excited states have 
been predicted by the extended Hubbard model including 
non-local interactions [TTHT5] . In particular a new kind 
of bound state between a holon and a doublon was re- 
cently introduced and named Hubbard exciton (HE) [TSl - 
[T9] . While the formation of HEs can be driven by a drop 
of Coulomb energy [20j [21], as in simple semiconductors, 
a kinetic energy loss could further stabilize the excitonic 
state in magnetic environments[22 -24 . Non- localized 
HEs have been extensively studied in the framework of 
high-temperature super-conductivity [25] [26] in relation 
with the proposed kinetic energy driven formation of the 
condensate [27H29] . 

In this letter we show that an extended HE pic- 
ture rationalize the optical properties of YVO3, a case- 
study for Mott-insulators TMOs. In this scenario, the 
optical transitions (see FIG. 1) observed at 2.4 eV re- 
flects the "single-particle" band (SP) whereas the one 
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at 1.8 eV is attributed to a kinetic energy related Hub- 
bard exciton [30]. Our pump-probe spectroscopic mea- 
surements in the 1.65-2.75 eV range revealed that the 
spectral weight (SW) is directly transferred between the 
two peaks, confirming the excitonic nature of the low- 
energy feature. Together with this, we measured both 
thermal and non-thermal effects and quantified the ki- 
netic energy contribution to the formation of the HE. 
Our study provides a new methodology, based on both 
static and time-domain spectroscopy, that can be used to 
unravel the complex nature of high-energy excitations in 
insulating TMOs and, more broadly, to study the kinetic 
energy-based mechanisms in strongly-correlated materi- 
als. 



The observation of multiple temperature-induced mag- 
netization reversals [31. accompanied by a series of 
structural [32 , magnetic [33 , and orbital transitions [34] 
turned the Mott insulator YVO3 into a case-study among 
TMOs. At room temperature YVO3 has a Pbnm or- 
thorhombic unit cell, where each V06 octahedron is 
tilted and distorted. The V 3+ ions have a 3d2 electronic 
configuration [35] so that two electrons occupy the t^g or- 
bitals. At 200 K a phase transition to a monoclinic phase 
(P611[5B] or P2i/c[32l [371 [38]) occurs with G-type or- 
bital order (G-OO). In this state, the d yz or d zx orbitals 
are alternately occupied both in the ab plane and along 
the c axis (FIG. 2 (e)). While still monoclinic and G-OO, 
at Tso=116K a C-type spin order emerges, character- 
ized by an antiferromagentic (AFM) spin configuration 
in the ab plane and ferromagentic (FM) along the c axis 
(FIG. 2 (f)). For further cooling below T^ ei)2 =77K a 
structural, orbital, and magnetic transition occurs: the 
system recovers the Pbnm crystalline structure while the 
orbital ordering switches to C-type (C-OO) and the spin 
order to G-type fG-SO)[39ti42]. The low-temperature 
phase will not be further discussed in the following (see 
additional materials). 
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FIG. 1: Reflectivity as a function of temperature of 
YVO3 obtained by ellipsometry data (a). The ratio 
between the spectral weight of the excitonic and single parti- 
cle band (b) increases rapidly while entering the G-00 phase 
and further increases at the spin ordering temperature. Both 
the SP and the HE spectral weights raise upon cooling (c). 
In (b) and (c) the dashed lines represent the transitions tem- 
peratures towards the orbital ordering (Too =200 K) and the 
additional spin ordering (Tso=H6 K). In (b) the ratio is nor- 
malized to the value at 300 K. Note that the reflectivity mea- 
surements in (a) are displaced for clarity from the measure- 
ment at 80 K. 



I. RESULTS 

All phase transitions can be identified by monitoring 
the optical properties of YVO3. The region of the optical 
spectra that is particularly sensitive to the orbital physics 
is the visible range. The two peaks characterizing the op- 
tical properties with energies of 1.8 eV (HE) and 2.4 eV 
(SP) are largely debatedgQl [42j OHSJ 02] and com- 
monly assigned to dfd^-djd 1 - transitions between two dif- 
ferent Vi and Vj sites[36] [43] . Multiplet calculations [36] 
account for the presence of the different peaks but fail 
in reproducing their temperature dependence: both the 
1.8 eV and 2.4 eV excitations gain spectral weight with 
decreasing temperature and approaching the spin order- 
ing transition at T5o=116K (FIG. 1 (c)), indicating that 
they correspond to the same high-spin state [30]. The S W 
gain with decreasing temperature is much stronger for 
the lower peak, which tentatively has been attributed to 
excitonic behavior. In fact from Hunds rules we expect 
the kinetic energy contribution to the formation of the 
Hubbard exciton to be more relevant in a FM-ordered 
phase[18 . 

In order to distinguishing the effects of temperature 
and orbital disorder we performed pump-probe reflectiv- 
ity measurements. 

The complex dielectric function e(uo) was measured by 
ellipsometry, for details see Ref.[30]. The static normal- 



incidence reflectivity R(lj) reported in FIG. 1 was cal- 
culated from s(lj). Broadband super-continuum probe 
experiments combined with an ultrafast optical pump at 
1.55 eV were performed on freshly polished ac oriented 
YVO3 samples mounted on the cold finger of a helium- 
flow cryostat. The reflectivity changes as a function of 
pump-probe delay ^(uj,t) = ^^^T^^ induced by 

80 fs pump pulses (E purnp =1.55eV, fluence<4mJ/cm 2 , 
at 40 KHz repetition rate and with polarization parallel 
to the a axis) were measured as a function of energy (for 
1.65 eV< Eprobe < 2.75 eV) and temperature. The lin- 
earity of the response was checked in all phases up to 
8 m J / cm 2 . 

The three phases are characterized by the different re- 
sponses summarized in FIG. 2 for probe polarization par- 
allel to the c axis (see online additional material for the 
intermediate temperatures). 

At room temperature, the pump-probe measurements 
are solely characterized by a very fast negative varia- 
tion of the reflectivity (decay time r ~ 0.5 ps) extending 
over the investigated spectral region (FIG. 2 (a)). On 
the contrary, the low-temperature phases are character- 
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FIG. 2: Broadband transient reflectivity spectra. 

Transient reflectivity as a function of wavelength (energy) 
and pump-probe delay for the three different phases: (a) dis- 
ordered, T=300K; (b) G-type orbital order, T=140K; (c) 
C-type spin order, T=80K. Gray lines report the transient 
reflectivity at the fixed energy of 2.33 eV as retrieved in stan- 
dard single-color pump-probe measurements. The oscillatory 
trend has been assigned to acoustic vibrations [46] and will be 
ignored. For each panel, the respective ordering patterns are 
sketched on the right hand side (d, e, f). 
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FIG. 3: From transient reflectivity to spectral weight. 

The transient reflectance is fitted using a variational approach 
on the static fit. (a), (b) and (c) depict the transient re- 
flectance for two characteristic times, typical of the "fast" 
and "slow" dynamics, at 300 K, 140 K, and 80 K, respectively, 
and the relative variational fit. In the inserts the static re- 
flectance for the different phases are reported 30 . (d), (e) 
and (f) show the time evolution of the SW variations in the 
different phases (see text). 



ized by a more composite optical response with fast and 
slow components, confirming previously reported single- 
color measurements [46 , 47 . The time-domain reflectivity 
measurements in the G-00 phase (116 K < T < 200 K) 
are characterized by a slow negative response centered at 
1.94 eV (FIG. 2 (b)), which gets more pronounced as the 
temperature is lowered (FIG. 2(c)). Finally, entering the 
S-00 phase (80 K < T < 116 K) a positive variation of 
the reflectivity appears at energies higher than 2.1eV. 



II. METHODS 

Our time-domain spectroscopic technique offers a di- 
rect view on the pump-induced changes of the reflectivity 
over a broad frequency range. This has significant advan- 
tages compared to a single-color pump-probe experiment. 



It allows, for instance, to determine the time-dependent 
spectral weight of the different features characterizing the 
optical response. The observation of a direct transfer of 
SW between the HE and SP band in YVO3 is the central 
experimental result of this study. 

In order to calculate the pump-induced evolution of 
the SW of the two bands from the transient reflectivity 
data, we proceed as follow. We consider the model used 
to fit the static ellipsometry, we calculate the equilib- 
rium reflectivity (Rq) and we fit the measured transient 
reflectance (AR exp (t) / R exp ) with a differential model 
(R(t) — Ro)/Ro, where R(t) is a model for the perturbed 
reflectivity obtained by variation of the parameters used 
to fit the equilibrium data. 

Among the 21 parameters used to fit the features in 
the static e(u) (six Gaussian and one Tauc-Lorentz os- 
cillators, see additional materials), the minimal set of 
parameters needed to obtain good differential fits at all 
temperatures and times contains only the following four: 
the amplitude, central frequency, and width of the os- 
cillator describing the HE and the amplitude of the SP 
one. FIG. 3 (a), (b), and (c) show typical fits obtained 
for the fast (blue) and slow response (red) in the different 
phases. The obtained values for the temporal evolution of 
the oscillator parameters are used to calculate the time- 
domain evolution of the spectral weight. The variation 
of the spectral weight for the two bands in time is cal- 
culated independently as the difference between the SW 
calculated for the model at a specific time (SW^^$ e p(t)) 
minus the static spectral weight (SWf}^ U g P ) (each spec- 
tral weight is calculated by numerical integration of the 
optical conductivity describing the band). The time de- 
pendence of ASW(t) = SW pumped (t) - SW static of HE 
and SP as a function of the pump-probe delay "t" is re- 
ported in FIG. 3 (d), (e), and (f) for three significative 
temperatures. The fast response (t < 3ps) in all phases 
is entirely described by a variation of the SP peak while 
the HE seems to be unaffected by photo-excitation in the 
first few picoseconds. This evidence can be rationalized 
as ground state depletion, confirming that only the higher 
energy optical transition is of single-particle origin. On 
the other hand the slow response (t > 10 ps), revealed in 
the orbit ally ordered phases, is related to spectral weight 
changes of both HE and SP bands. 



III. DISCUSSION 

The time-domain response cannot be accounted for by 
photo-induced heating. The spectral weight of both oscil- 
lators increases upon cooling (FIG. 1) and it is therefore 
expected that a transient laser-induced heating would re- 
sult in a decrease of the spectral weight of both SP and 
HE. However, the time-domain measurements (FIG. 2) 
reveal that only the SW of the HE peak decreases, while 
the SP band shows the opposite behavior, which rules out 
a simple heating effect. The non-thermal SW gain of the 
high energy oscillator (FIG. 4) lasts up to ~ 400 ps and 
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FIG. 4: Long timescale pump-probe measure at 80 K. 

a) Time-domain reflectivity data at 80 K. Note that the fast 
response visible in FIG. 2 is absent beacuse of the coarse tem- 
poral step used, b) Transient reflectance at t = 40 ps (red 
curve) and 900 ps (black curve). The expected thermal re- 
sponse (R static (85K)-R static (80K))/R static (80K) is shown 
for comparison (green curve). 



only at longer times the measurements indicate a SW 
loss for both oscillators. A comparison between green 
and black curves in FIG. 4b reveals that 1 ns after pho- 
toexcitation the degrees of freedom have not reached the 
thermal equilibrium. 

Following those considerations we analyze the transient 
spectral weight measured at t > 40 picoseconds (well be- 
yond the electronic relaxation) as resulting from a ther- 
mal contribution and a non-thermal one. The thermal 
contribution to the SW variation of the HE and SP 
peaks (SWg% tic (T + AT) and SWgp tic (T + AT)) can 
be calculated by interpolation of the temperature be- 
haviour of the static measurements at T + AT, where 
AT is the photo-induced heating obtained from thermo- 
dynamic considerations (details in additional material). 
The SW variations of non-thermal origin can therefore 
be calculated subtracting the thermal contribution to the 
experimental values: 

ASW™£- thermal SWff J £ iped (50 ps) - SWhe Uc (T + AT) 

ASW n s ™- thermal = SW^ mped (50ps)-SW^ tic (T+AT), 
for HE and SP respectively, where SWf umped (50 ps) (i = 
HE, SP) are the means of the measured photo-excited 
spectral weights at t=50±30 ps after the pump arrival. 

The non-thermal components of the pump-driven spec- 
tral weight variations at different equilibrium temper- 



atures are reported in FIG. 5. It is evident that the 
non-thermal contributions consist of a direct exchange 
of spectral weight between the HE and SP: this result 
prove that a the two lowest lying optical excitations be- 
long to the same multiplet and that the 1.8 eV transition 
is of excitonic nature [30]. Moreover, the photo-induced 
increase of SW transfer from the HE to the SP in the 
C-SO phase (FIG. 5 for T < 120 K) highlighs the kinetic 
energy contribution to the formation of the HE. 

We argue that the spectral weight loss of the HE in 
favor of the SP is driven by pump-induced spin-disorder 
(FIG. 6). Two adjacent and c-oriented orbital-chains of 
YVO3 in the G-OO/C-SO phase are sketched in the first 
row of FIG. 6. The spin order is FM along the c axis 
(J c < 0) while it is anti- ferromagnetic in the orthogo- 
nal direction (J a > 0). The photo-exctitation (a axis 
polarized) transfers charges between chains, leaving one 
excess electron on a site and one vacancy on the near- 
est as depicted in FIG. 6 (t=0). The magnetic cou- 
pling between an excited and a non-excited site along 
c changes as a consequence of such a charge redistribu- 
tion (J*). This dramatic perturbation of the spin cou- 
pling can be grasped by the following considerations. The 
Goodenough-Kanamori rules [48 , stating that the super- 
exchange coupling J between half-filled and empty or- 
bitals is negative, is consistent with the FM order ob- 
served along the c axis. At simplest, the pump driven 
changes in the orbital occupation leads to half filling, in- 
ducing a switch from FM to AFM exchange coupling. 
The electronic subsystem relaxes within a few picosec- 
onds (ti) while the spins remain locally perturbed (tilted 
blue arrows in FIG. 6 for t=ti). We argue that, being 
the orbital angular momentum completely quenched, the 
coupling mechanism allowing the magnetic subsystem to 
transfer energy to other degrees of freedom is very weak. 
Therefore, the local perturbation of the spins diffuses (t2) 
leading to spin disorder. This partially suppress the ki- 
netic energy gain of the HE and lead, subsequently, to a 
non-thermal spectral weight transfer to the SP peak. A 
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FIG. 5: Non-thermal spectral weight changes of the 
HE and SP peaks. The excitonic nature of the low-energy 
transition is revealed by the direct SW exchange between the 
two excitations (see text). The error bars are estimated from 
the fitting procedure. 
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fully thermalized state is achieved only on much longer 
timescales (> Ins) through spin-lattice coupling [47]. 






FIG. 6: Non-thermal spin disorder. YVO3 is in the 
G-OO/C-SO phase; t indicates the pump-probe delay (t<0 
represents the unperturbed state). The photo-excited holons 
and doublons perturb locally the magnetic coupling along the 
FM chain J* (J c '*) (t=0, central site). The photo-excited 
electrons relax within few picoseconds (t~ti) leaving a local 
perturbation on the spin system. The spin disorder diffuses 
along the FM chains on a longer timescale (t=t2) resulting in 
SW transfer between the HE and SP peaks (see text). 

In conclusion, making use of a novel time-domain 
methodology, we show that the optical properties of 
YVO3 can be rationalized introducing a bound state be- 
tween a doublon and a holon, named Hubbard exciton 
(HE). In this model the two oscillators at 1.8 eV (HE) 
and 2.4 eV (SP) are described as a single band originat- 
ing from the same inter- Vanadium transition. Moreover, 
we propose that the time-domain measurements of the 
spectral weight exchanged between the HE and the SP 
band can be used to study the kinetic energy gain of the 
excitonic feature associated to spin order. Further stud- 
ies are called to verify these new intriguing scenarios. 

The authors are grateful to M. Malvestuto, J. van den 
Brink, L. Hozoi and K. Wohlfeld for useful discussion. We 
acknowledge G. R. Blake for reviewing the manuscript. 
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Supplementary information 

IV. STATIC FITS 

We fitted 62 by a sum of six Gaussian peaks and a 
Tauc-Lorentz oscillator for the transition at lower energy 
(HE), ei by the Kramers-Kronig consistent functional 
form. With those line-shapes, that have been justified 
and extensively used to address amorphous and locally- 
disordered materials [1, 2 , we obtain static fits of very 
good quality (FIG. 7). 

These are the dielectric expressions used: 

e 2 (oo) = L a (oo) + Lb(u) + . . . + L g (cj), 

where L a (oo) = ^^r^p^^j if 00 > gap and 
L a (uu) = elsewhere. L c , L e , Lf, L g are analogous 

of Lb(uj) = b(e 7 fc — e 7fc ); a,b,c are the 

amplitudes, 00 i and 7^ the central frequencies and widths 
for i — a, 6, c. . . 

ei(cj) = 1 + D a (oo) + D b (oo) + . . . + A,M, 
where D a (oo) has the Jellison's form[3] and the other 
terms are proper sums of Dawson's functions [4]. By si- 
multaneous fitting of both the real and imaginary part 
of the dielectric constant we were able to obtain a series 
of continuously-varying parameters. In the c-axis fitting 
the Gaussian centered at 33000 cm -1 fades below 140 K: 
we left this oscillator free to vary as attempts to block 
it produces a much worse fit. For the same reason we 
bind the central frequency of another oscillator at 23000 
cm -1 . These assumptions are justified by the overall 
good fits of the static optical properties. It should be 
noted that, this approximations, which have been made 
to have a good description of the static optical properties 
outside of the interest range for the time domain data, 
do not affect the outcome of the time domains fits. This 
was verified by performing time domain fits with differ- 
ent sets of parameters. The spectral weights have been 
computed as f^ 3000 00 62(00) doo, with 00 in cm -1 and €2(00) 
composed only of the HE peak or the SP peak. 



V. TIME-RESOLVED MEASUREMENTS 

We measured the transient reflectivity AR(t)/R as a 
function of temperature and pump-probe delay in the 
450-750 nm wavelength-region after excitation with 4 
mJ/cm 2 of 775 nm ultra-short Ti:Sa laser pulses at 40 
KHz. The linearity of the response was checked in 
all phases. Time-domain measurements were performed 
with the pump parallel to the a-axis (P||a) and probe 
parallel to c (p||c) for the intervals -2-1-4 ps and 4-1-80 ps. 
A set of representative measurements is plotted on FIG. 
8. Similar measurements performed with pump parallel 
to the c-axis (P||c) and probe parallel to the a-axis (p||a) 
shows only a temperature-independent fast decay time 
(FIG. 9). Long timescale measurements for P||a and p||c 
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FIG. 7: Imaginary (top) and real (bottom) part of the di- 
electric constant of YVO3 measured by ellipsometry[30 along 
the c axis at 80 K and 300 K. The vertical dashed lines mark 
the positions of the central frequencies of the oscillators used 
to perform the static fitting procedure (see text). 

are shown in FIG. 10. The data shown are limited to 
temperatures as low as 80 K, because the crystal tend 
to break at the low temperature phase transition often 
leading to a loss of thermal contact [30] and strong aver- 
age heating. 

VI. ANALYSIS OF THE SPECTRAL WEIGHT 

Starting from the fits to the static optical properties 
and our pump-probe data we extract the time-domain 
spectral weight variation ASW(t) of the HE and SP ex- 
citations as follows. 

We fit a differential model (R(t) — Ro)/Ro to the mea- 
sured transient reflectance AR exp (t) / R expi where Rq and 
R(t) are the models describing, respectively, the equilib- 
rium reflectivity and the pump-perturbed one as a func- 
tion of pump-probe delay t. The values of the oscillators 
parameters obtained by this fitting procedure at different 
times are used to calculate the evolution of the spectral 
weight: 

ASW(t) = sw pumped (t) - SW static (T), 
where SW pumped (t) is calculated from the fitting pa- 
rameters of R(t) and SW staUc (T) is obtained from the 
static reflectivity at temperature T (both SW are calcu- 
lated by the numerical intergration described earlier). 

A. Non-thermal contribution 

At any fixed temperature T, the non- thermal contribu- 
tion to the variations of the SW of HE and SP can be cal- 
culated from static optical properties, the time-resolved 
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FIG. 8: Relative variation of the reflectivity in the -2 to 80 
ps range, with P||a and p||c. 



data and the laser pump energy, as follows: 

/\SW non ~ thermal (t) = SW pumped (t) - SW static (T + AT(t)), 
where sW pumped (t) is the photo-excited SW and 
SW static (T + AT(t)) is obtained by interpolation at 
T + AT(t) of the static model. AT(t) is the pump- 
induced heating calculated from a two-temperature 
model (2TM)EH3 for the lattice (L) and spin (S) de- 
grees of freedom: 

C L ^t = - 1 {T L -T s )+pP efJ {t) 
Cs^=-l{Ts-T L ) + (l-p)P eff {t) 
where Cl and Cs are the heat capacities [32] of the 
two subsystems, 7 is the magnetoelastic coupling and p 
a phenomenological constant. In this model we assume 
that the pump pulse P(t) photo-excites carriers from the 
lower Hubbard band (LHB) to the upper Hubbard band 
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FIG. 9: Measurements with P||c and p||a. 
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FIG. 10: Relative variation of the reflectivity in the -40 to 
1000 ps range, with P||a and p||c. 



(UHB). As the quasi-particles relaxes, they act as an ef- 
fective pump P e ff(t) for the lattice and spin degrees. In 
this picture, p represent the phenomelogical coupling of 
the electronic subsystem to the other two. The behaviour 
of T L and T s is reported in FIG. 11 for T=80K. 

The validity of this model is confirmed by comparison, 
at any temperature, with the expected thermodynamic 
steady-state temperature increase AT(T). It is straight- 
forward to write 

AT(T\ — Qabs-N A -V ^ 150 

L ^ ± K- 1 ) — S-d-u-C mo i ~ C mo i[J/(mol-K)]i 

where Q a bs is the pump energy absorbed by the sam- 
ple, Na the Avogadro's number, V the elementary cell 
volume, S the sample's surface irradiated, d the pump's 
penetration depth, u the number of chemical units in a 
cell and C mo i the temperature-dependent total heat ca- 
pacity. 

There is a good agreement between the temperature 



8 



„ 30 










s — ■ — -^^L^" 1 







Spin 






— Lattice 








AT(T=80K,t^co) 





3456 23456 23456 

10 100 1000 

Delay [ps] 



FIG. 11: Two temperature model for T=80K. A gaussian 
P e ff(t) with FWHM = 3ps is turned on at t = 10 ps. T s (t) 
is reported in red, Tl (t) in green while the converging straight 
line AT(80K,t oo) is dotted. 



dependence of the temperature increases for the two 
models, as shown in FIG. 12 (red and black curves). At 
this point, the 2TM permits to obtain the temporal de- 
pendence of the lattice temperature and allows for the 
calculation of the non-thermal component. The blue dots 
in FIG. 12 represent the temperature variations at pump- 
probe delay t = 50 ps used to obtain the non-thermal 
contribution to the variation of the spectral weight (FIG. 
5). 
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FIG. 12: Comparison of the two thermodynamic models at 
equilibrium or long pump-probe delay. The blue dots repre- 
sent the calculated lattice temperature used to calculate the 
non-thermal contribution to the SW variations at different 
temperatures. 
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